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Abstract 

We consider various 2D lattice equations and their integrability, from the point of 
view of 3D consistency, Lax pairs and Backhmd transformations. We show that these 
concepts, which are associated with integrabihty, are not strictly equivalent. In the 
course of our analysis, we introduce a number of black and white lattice models, as 
well as variants of the functional Yang-Baxter equation. 



1 Introduction 

Recent progress in the description of integrable partial difference equations is to a great 
part due to the consistency approach [H El [3l H] , in particular in the form of 3 dimensional 
Consistency- Around-a-Cube (CAC). One of the highlights of this approach is the immediate 
existence of a Lax pair and Backlund transforms (BT), which can be directly constructed 
from the "side-equations" of the cube |21 [S]. One can say, in effect, that the set of side 
equations yields the BT and the Lax pair. One situation where side equations have been 
used effectively is in constructing soliton solutions to the lattice equations [HI [3, IH]. 

Originally it was assumed that on all faces of the cube the equations were the same 
in form, depending on the relevant corner variables (one component at each corner) and 
spectral parameters. Recently, 3D consistent sets have appeared, with different equations 
on the faces O [101 HH [121 [13] . In this more general context it is of interest to take a closer 
look at the Lax pairs, Backlund transforms and consistency, and investigate what they are 
good for. It should be noted that in the context of partial differential equations the existence 
of trivial Lax pairs is well known [1^ and similar examples have also been noted for some 
discrete equations (see [I^, Chapt. 6). 
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We will first recall how consistency around the cube, existence of Lax pairs and Backlund 
transformations are intimately related for lattice maps on the square lattice given by multi- 
affine relations (sections |2] and [s]). These considerations apply to the elementary cells, and 
are local. In section |4| we describe specific examples having the CAC property. We detail in 
particular the already known explicit forms of the equations and Lax pairs which we use in 
the rest of this paper. In particular, we find that in some cases the zero curvature condition 
(ZCC) yields two different equations that can be used to define rational evolution in the 
lattice. We then address the global problem of defining equations over the whole lattice, 
with the guideline given by the 3 dimensional structure coming from the CAC construction 
(section [5]), and check integrability with the calculation of algebraic entropy. In section [6| 
we push the use of 2 x 2-matrix Lax pairs to its limits, by constructing discrete systems over 
a larger sub-lattice of the original lattice. This brings to the light interesting (integrable) 
structures related to a generalized form of the functional Yang-Baxter equations pjj . 



2 3D consistency, a reminder 



The starting point is a regular 2D square lattice, with vertices labeled by integers n, m. Func- 
tions Xn^m are associated to the vertices, and they are subject to a constraint at all elementary 
cells. This constraint is expressed by an equation Q{xn,m,Xn+i,m,,Xn,m+i,Xn+i,m+i) = 0, as- 
sumed to be multi-affine in the four vertex variables. It should depend on all 4 vertex 
variables, and it should not factorize. It may also depend on some parameters. Sometimes 
the parameters can be associated to specific directions of the lattice, in which case they ap- 
pear as "spectral parameters" . To ease the notation, one usually denotes the running value 



Xn. 



by Xr, 



X. 



in 3D setting, xi 



and for neighboring values one only indicates the shifts: 

= Xn+l,m,ky ^113 = ^n+2,m,fc+l CtC. 



X. 



n,m+l 



X2, or 



For multidimensional consistency one needs to build a cube on top of a square and give 
equations on all six faces of the cube, see Figure [1} the bottom equation being the original 
one. 



Xl3 



X3 - 






X , 
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' 3^12 



X2^ 



X2 



Figure 1: The consistency cube 



Supposing x,xi,X2,x^ are given, one can compute Xi2 using the bottom equation, xi^, 
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using the left side equation, and X23 using the back side equation. We may then calculate the 
final value X123 in three different ways, using the front-, right- and top-equations, respectively. 
The equations are said to be Consistent-Around-a-Cube (CAC) if the three ways yield the 
same value for 2:123. 

Note that it is also possible to check CAC with other initial values that allow full evolu- 
tion, for example, X3,x,Xi,Xi2- 



3 Lax ~ Backlund 



The Lax/BT approaches differ from CAC in that only the side equations are used as input, 
while the top and bottom equations are supposed to be derived from them. 



3.1 Construction of the Lax pair 

The Lax pair for the bottom equation is constructed from the side equations by isolating 
the variables xs, X13, X23, a;i23, and writing the equations as: 

left: xi3X3Ci(xi,x) +xi3C2(xi,x) +a;3C3(xi,x) + C4(xi,x) = 0, (la) 

right a;i23a;23Cl(Xi2,X2) +Xi23C2(Xi2,X2) +X23C3(Xi2,X2) +C4(X12,X2) =0, (lb) 

back X23X3bi{x2,x) + X23b2{x2,x) + X3bs{x2,x) + b4{x2,x) = 0, (Ic) 

front a;i23a;i3&i(a;i2,a;i) + X123 &2(a;i2, ^i) + Xi3b3{xi2,xi) + b4^{xi2,xi) = 0, . (Id) 
The coefficient functions bi,bi,Ci, Ci will be affine linear in their arguments. 
We now introduce the homogeneous coordinates /, g for X3 and its shifts by 

X3 = f/9, X23 = 12/92, Xr3 = fl/9U a;i23 = /l2/fi'l2- 

This amounts to considering 0:3 as belonging to the projective space CP^. We denote by ip 
the pair 




which is defined up to a global factor. Equations ([T]) can then be written as 

iIji ~ Lip, {iIj2)i ^ L1IJ2, il)2 ^ Mil), (^1)2 ^ M{il)i), 

where 

L{x,,x) ^ f-^3(xi,x) -c,{x,,x)\ ^ f-Y^^'f -Mx2,xJ\ 

\Ci[Xi,x) C2[Xi,x) J \bi[X2,x) b2{X2,X) J 

and similarly for the bar-quantities. Since we are working in CP^ all equalities are projective 
equalities, and as a reminder of this we use the symbol ~ to indicate that two matrices are 
equivalent if their entries are proportional, in other words, L and M belong to PGl{2, C). 
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The Lax matrices provide parallel transport of i/j along the bonds of the lattice. The zero 
curvature condition means that the parallel transport along any closed path on the lattice 
is trivial. It is necessary and sufficient to ensure that taking from position (0, 0) to (1, 1) 
via the two routes (0, 0) — > (0, 1) — > (1, 1) and (0, 0) — )■ (1, 0) — )■ (1, 1) gives the same result: 

(?/'i)2 ~ (V^2)i, i.e., L(xi2,X2)M(x2,x) ~ M(xi2,Xi)L(xi,x). (2) 

This matrix relation yields three scalar equationsj^ These equations are written in terms of 
the variables x, Xi, 0:2, Xi2, and in order to satisfy the ZCC we must impose a constraint on 
these variables. In standard cases, this constraint will just be the bottom equation. Indeed, 
in integrable cases the three equations in ^ have a common factor. It may also happen that 
the ZCC is satisfied automatically, if the side equations are simple enough, or sometimes the 
common factor may factorize, yielding multiple choices. We will examine specific examples 
of this phenomenon below. 



3.2 Direct approach: Backlund transformation 

A Backlund transformation (BT) is generated by a set of equations S{X, Y) = on two sets 
of functions X and Y, such that eliminating one set of functions say X (resp. Y) one gets 
an equations E{Y) = (resp. E'{X) = 0) on the other set (for details, see Proposition 4.1 
in p.) 

For a 3D consistent set of equations the side equations provide a BT between the top 
and bottom equations. In order to derive the bottom equation from the side equations one 
proceeds as follows: 



1. Solve Xi3 from the left equation (la) 



2. Solve X23 from the back equation (Ic). 

3. Solve X123 from the front equation (Id). 
4. 



After this the right equation (lb) can be written as 0:37^2 + x^TZi + TZq = 0, where TZi 
are polynomials in x, xi, 0:2, a;i2. The greatest common divisor (CCD) of TZi (or one of 
its factors if it factorizes) will be the bottom equation. 



One can find the following explicit forms for the TZf 
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^One could do the computations in some particular representative of the equivalence class, e.g., by 
requiring the matrices to be uni-modular. However this is not necessary, and may in fact be cumbersome, if 
it introduces square roots. 
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That the Lax and the BT approaches yield the same equations can be seen as follows: 
Let us denote V := ML, U := LM, then the matrix elements of V, U are related to TZi as 
follows (here subscripts indicate the matrix element): 



U21 


V21 


7^2 


U22 


V22 


7^o 


Un 


Vu ~ 
U21VU 


VnUn' 
U12V21 


Uu 
U22 
Un 


V12 
V22 


V12U12 




UnVn 


UnVn 


Vn 


UnVn 



Thus f/ ~ V iff 7^^ = 0. 

Of course one can equally well use the side equations to solve for the variables xi, X2, X12 
from the left-, back-, and front-equations, after which the right-equation will be a polynomial 
in X with coefficient depending on X3, Xis, X23, 2:123, with their GCD yielding the top equation. 



4 Examples 

4.1 Linear side equations 

As a first example we consider the case where all the side equations are linear, i.e., the left 
side equation is: 

xi'i - xi - + X = Q (7) 

and we have the same equation with suitable subscript changes on the other vertical sides. 
One then finds that neither the Lax nor the BT approach yields anything. The conditions 
are satisfied without reference to the bottom equation. 

What about GAG? It will involve both the bottom and the top equations, so it should 
give some conditions. Indeed, if one tries GAG with ([t]) and a completely general multi- 
affine bottom and top equations, related by a shift (but with same parameters) one finds a 
consistent set with the bottom equation 

a{x - Xi){x2 - X12) + h{x - X2){xi - X12) + c{x - Xi- X2 + x^) + d = 0. (8) 

This is a combination of the lattice modified KdV equation (ImKdV) (aka Ql in the ABS 
list [1]), and the linear equation. 

A similar analysis can be done starting from the linearizable side equation 

X13X = X1X3. (9) 

The Lax matrix for this system is diagonal and the ZGG is automatically satisfied. From 
GAG analysis we find that ^ is compatible with the six parameter family of homogeneous 
bottom/top equations of degree 2 

ai xx\ + a2 XX2 + 03 XX12 + a4 xi X2 + 05 xi x\2 + X2 X12 = 0. (10) 



We will examine the integrability of these equations in section (5.1) 
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4.2 HI 



The lattice potential KdV (IpKdV), which describes the permutability property of continuous 
KdV, is a paradigm of integrable lattice equations (aka HI in the ABS list [4j). For this 



system everything works well. [We will return to this model in a new context in section 6.2 
The model is given by 

B.1 := {xi - X2){x - X12) - p + q = 0. (11) 

After imposing this equation on the sides of the CAC cube, with suitable variable changes, 
i.e., 



left 
right 
back 
front 



(x - Xi3)(xi - X3) - - r) = 0, 
(x2 - a;i23)(a;i2 - 2:23) - (p-r) =0, 
{x - X23){x2 - X3) - (g - r) = 0, 
(xi - Xi23)(a;i2 - 3:13) - (g - r) = 0, 



one easily finds 



= (a;i2-x)-Hl, 
TZi = [(x - Xi2)(xi + X2) - a - /3 + 27] • HI, 
TZo = [xiX2(xi2 - x) + (/3 - 7)xi + (a - 7)^2] ■ HI, 



(12a) 
(12b) 
(12c) 
(12d) 



(13a) 
(13b) 
(13c) 



with HI, as given in (11), as the GCD. Similarly by working on the bottom variables 



x,xi,X2,xi2 one obtains the 3-shifted HI as the top equation. 
The Lax matrices in this case are 



L(xi,x) = 
and one easily finds that 



X p — r — xx\ 
1 — Xi 



M(X2,X) 



X q — r — XX2 
1 -X2 



(14) 



M(xi2, xi) L(xi, x) - L(xi2, X2) M(x2, x) = HI X 



1 -(X1+X2) 
-1 



4.3 Hlgi a deformed version of HI 



This is an asymmetric deformation of ( [l2| ) [9l |TT] 

left: (x — Xi3)(xi — X3) — (p — r)(l + exxis) = 0, (15a) 

right: (x2 - Xi23)(a;i2 - 3:23) - (p - r)(l + 6X12X23) = 0, (15b) 

back: {x - X23)ix2 - X3) - {q - r){l + exx23) = 0, (15c) 

front: (xi - X123) (2:12 - ^13) - (g - r)(l + 6X13X12) = 0. (15d) 

The BT or Lax computations give as GCD the bottom equation 

(x - Xi2)(xi - X2) - (p - g)(l + exxi2) = (15e) 
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and similarly for the top equation we get 

{x3 - Xi23)(a;i3 - 3:23) - (p - + £3:13X23) = 0. (15f) 

If we draw a line connecting the corners appearing in the deformation term as a product, 
the lines form a tetrahedron inside the cube, see Figure |2} 




Figure 2: The cube of Hie 

The special thing about these equations is that the parallel sides are not obtained by 
shifts, but the shift must be accompanied by a 90° rotation, as described in Figure |2j The 
fact that the parallel sides are not identical implies that the Lax matrices associated to 
parallel bonds come in two different forms. Letting 

L{xi,x)=\ ^ ) ' M{x2,x) = L\p^g^^^^^^, (16) 

and 

L'ixi,x) = ' M'{x2,x) = (17) 

we have 



M'{xi2,xi) L{xi,x) - L'{xi2, X2) M{x2, x) = 



IM X 



1 - 6X1X2 -(xi + X2) 
— exiX2 —1 + exiX2 



Since we have two different parallel sides we can have two different bottom equations, one 
derived from M'(xi2, xi) L(xi, x) — L'(xi2, X2) M(x2, x) as above, the other possible choice 
being 

'1 -(X1 + X2)' 



M(xi2,xi) L'(xi,x) - L(xi2,X2) M'(x2,x) = (15e) x 







-1 



The totality of equations on this cube, obtained from the previous one with L ^ L' etc 
actually corresponds to an inversion of the cube by 



X O X123, Xi ^ X23, X2 O Xi3, X3 ^ X12, 
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which yields 



left 


(x — 


Xl3){xi — X3) — 


{p 


— r)(l + exiXs) = 0, 




right 


- 


- Xi23)(Xi2 — X23 


)- 


- r)(l + 6X2X123) 


= 0, 


back 


{x — 


3^23) (3:2 - 3^3) - 


(q 


- r)(l + 6X2X3) = 0, 




front 


(xi - 


- a;i23)(3^12 — 3^13 


)- 


(g - r)(l + exiXi23) 


= 0, 


bottom 


(x — 


Xi2)(Xi - X2) - 


{p 


- g)(l + 6X1X2) = 0, 




top 


{X3 - 


- 3:123) (a;i3 - a;23 


)- 


{P - 9)(1 + £2:3X123) 


= 0. 



(18a) 
(18b) 
(18c) 
(18d) 
(18e) 
(18f) 



These are also consistent. We thus have two kinds of cubes, one defined with (15), and one 
with (18). These cubes can be glued together: note that the 2-shift of (18a) is equal to 



(15b). In fact under the central reflection the diagonal lines in Figure 2 would change to the 



other diagonals. 

We will return later to the question of gluing the cubes to fill the space. 



4.4 Flipped HI, 



This model was proposed in [13] (see Sec. 3.1, Case (e, 0, 0, e)). In that paper the model was 
given in a cube with flipped coordinates X2 -H- X23, xi -H- X13. In the coordinates of Figure [T] 
the side equations of this model are given by 



left 
right 
back 
front 



(x - Xi)(xi3 - X3) - {p- r)(l + exxi) = 0, 

(X2 - Xi2)(Xi23 - 2:23) - (p - r)(l + eX2Xi2) 

(x - X2)(x23 - 3:3) - (g - r)(l + exx2) = 0, 

(Xi - Xi2)(xi23 - 3:13) - (g - r)(l + 6X1X12) 



(19a) 
(19b) 
(19c) 
(19d) 



BT derivation of bottom and top-equations. Now computing the values of X13, X23, X123 
from left-, back- and right-equations, respectively, gives for the front equation an expression 
that does not contain X3 at all. This expression factorizes into two factors and thus we could 
have two different bottom-equations: 

bottoml exxiX2Xi2(l/x — 1/xi — I/X2 + I/X12) — x -|- xi -|- X2 — X12 = 0, (20a) 
bottom2 p{x — X2)(xi — X12) + g(x — Xi)(xi2 — X2) + r(x — Xi2)(x2 — Xi) = 0. (20b) 

Similarly, working with the values at the bottom square we get two candidates for the top 
equations: 

topi X3 - Xi3 - X23 + 3^123 = 0, (21a) 

top2 p(x3 - X23)(xi3 - X123) + q{x3 - a;i3)(a:i23 - X23) + r(x3 - Xi23)(a;23 - 3:13) 

+e{p — q){q — r){r — p) = (21b) 

Notice that in down2 and top2, there is an explicit dependence on r even though these 
are 2D equations. In fact, one should now consider r as a global parameter, although it was 
associated with the third dimension in the above derivation. 
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Performing the usual CAC computations with the given sides and different top and 
bottom equations reveals that the set of equations is consistent in two cases: with the pair 
(bottoml,topl), or with (bottom2,top2), which was given in |9]. In other words, given the 



side equations ( 19 ) there are two consistent ways of completing the cube. Continuing further, 
we can also interpret the cube to provide a BT between the left and right equations. Indeed, 
if we do the above BT construction on the corners of the left side equation the result is 
identically zero with (bottoml,topl), while (bottom2,top2) produces the right-equation. 



The Lax matrices. The standard procedure gives 

T, \ A , A A((7,x, X2)\ w \ (a - r)(l + exy) 
L(a;i,x)=L ^^'^''M, M(x2,a;)=L ^^'^''M, A(a, a;, y) := ^ ^ ^ 

^22) 

Since the matrices are upper triangular the ZCC implies 

S := A(p, X2, xvi) + A(g, x, x-^) - A(g, xx.xvi) - A(p, x, xi) = 0. (23) 

Remarkably enough, the above sum factorize as 

S = bottoml ■ bottom2 /[(a; — x\)[x — X2)[x\ — X\2){x2 — x\2)\. 

Note that we can write ( [23| also in the form 

S = (T - l)A(p, x, xi) - {S - l)A(g, x, X2) = 0, 

where T is a shift in m and S a shift in n. This of course is in the form of a conservation 
law. 



5 Filling the space with consistent cubes 



So far we have only considered a single cube and its CAC/Lax/BT. But as the name indicates, 
lattice equations should be defined over the whole lattice. This brings further complications, 
for example with one cube we could freely do different Mobius transformations in each corner 
of a cube, but when the cube is part of lattice such seemingly innocuous actions will affect 
neighboring cubes as well and can destroy the lattice structure. 

The rule is simple: cover the two dimensional lattice with consistent cubes, with the 
condition that adjacent vertical faces coincide exactly, that is to say their four corner values 
satisfy the same equation. This is expected to produce integrable lattice equations. 

We will follow this guideline for various models mentioned before, and systematically 
check integrability of the lattice equations so obtained, by calculating their algebraic en- 
tropy [I71[T8]: the vanishing of the entropy is a yes/no test which gives a clear cut separation 
between integrable and non integrable cases. 

We may briefly recall how to calculate the entropy. The local equation determines an 
evolution, starting from initial conditions given for example on a diagonal staircase (lattice 
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points of coordinates {m,n) with m + n = or 1). The solution is then calculated on 
diagonals moving away from the diagonal of initial conditions, explicitly in terms of these 
initial conditions. The algebraic entropy is defined as the rate of growth of the degrees on 
these diagonals. Exponential growth is generic, and polynomial growth is characteristic of 
integrability, while linear growth is associated with linearizable equations p!7l [T8| [19] . 

The exact shape of the diagonal line on which the initial values are given is not important. 
If one modifies this shape locally, the sequence of degrees will change, but not its asymptotic 
rate of growth. This should be kept in mind for some of the models studied below. For 
example one could very well change the initial diagonal (steps of height 1 and width 1) to 
a diagonal with bigger steps, for example with height 2 and width 2, as this will not affect 
the calculation of the entropy. 



5.1 Equations consistent with linear sides 



We can check the integrability of the two quad equations given in section (4.1). The first 
one, i.e., equation ([8|, leads to the sequence of degrees: 



{dn} = 1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 67, 79,... 

that is to say c/„ = 1 + (n^ + n)/2. This quadratic growth confirms integrability. 

Remark: Since it is integrable, it should also be consistent with a nontrivial set of side- 
equations, such that one can produce it via Lax/BT computations, but we will leave this 
open. 



For the second equation of section (4.1), that is to say the general homogeneous relation 



(10) of degree 2, the result is different. We get the following sequence of degrees: 

{dn} = 1, 2, 4, 9, 21, 50, 120, 289, 697... 
which may be fitted with the generating function 

,w = x:-'..^" (XV. -1) 

indicating a non vanishing entropy e = log(l + \/2), showing non-integrability. Thus even 
CAC is not sharp in this case. It was actually shown in [20] that the simple additional 
condition renders ^ integrable. 



5.2 Hie: Checkerboard lattice 

We have already noted that in the Hie model there are two different cubes related by 
inversion, and that these cubes can be glued together in a unique way, providing a black- 
white lattice. This problem has been discussed in detail in [TT]. 
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Since the gluing process is completely fixed and each cube has the CAC property it 
is expected that the composite lattice is integrable. We have calculated the sequence of 
degrees {dn} for the evolution defined by this model, with initial conditions on a diagonal 
as prescribed above. The outcome is the sequence 

{dn} = l, 5, 13, 25, 41, 61, 85, 113, 145, 181, 221,... 

fitted by the generating function 

The above sequence has quadratic growth. Thus the entropy is vanishing and the model is 
integrable. 



5.3 Flipped Hie: More black and white lattices 



In the flipped Hie case (Sec. 4.4), the side equations are the same, allowing simple gluing 



together of the cubes, but, as we saw in section ( 4.4[ ), the side-equations are somewhat weak 



and allow two different compatible pairs of bottom/top equations. We may then cover the 
two dimensional lattice with consistent cubes, assigning either equation bottoml (topi), 
which we call white, or equation bottom2(top2) which we call black to each elementary cell. 
This can be done in an arbitrary way if one just insists on having a compatible 3D structure 
of cubes over the 2D lattice. It is then natural to ask which of the configurations obtained 
in this way are integrable. 

Let us consider periodic distributions. The lattice is divided into rectangular groups of 
cells of width h and height v. Within such a rectangle, a fixed assignment is made, and the 
pattern is repeated periodically in both directions. (A pattern with v = 1 and h = 1 gives a 
uni-colored assignment.) 

Consider for example {h,v) = (2,2). There are a priori 2^ possible patterns of that 
size, but only 3 inequivalent ones which cannot be reduced to configurations having smaller 
periods (see Figure [3]). The naming convention is to list the colors starting from the lower 
left corner onwards, denoting bottoml/topl with 0, alias white, bottom2/top2 with 1, alias 
black. The equivalence of patterns comes from the fact that we have to look at the lattice 
globally. It is easy to see, for example, that in the case {h, v) = (2, 2) we have the equivalences 
[0100] ~ [0010] ~ [1000] ~ [0001], [1011] ~ [1101] ~ [0111] ~ [1110], and [1010] - [0101] 
(checkerboard lattice). Moreover [0000] and [1111] have periods (1,1), [0101] and [1010], 
[0011], [1100] have periods (2,1) and (1,2). 

Claim: Some of the distributions are integrable, and some are not. Although the pattern is 
3D consistent, the Lax pair is weak and cannot precisely fix the bottom and top equations. 



1x1 patterns (unicolor distributions). Both unicolor distributions {h,v) = (1, 1) have 
vanishing entropy. The purely white one is linear. The purely black one is non-trivially 



11 



5 E 9 

[0110] [1011] [0100] 



Figure 3: The three inequivalent (2,2) patterns. 



integrable, showing quadratic growth of the sequence of degrees 
{dn} = 1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 56, 67, 79, 92, 106, 121, 137, 154, 172... (25) 

1x2 patterns. Both 1x2 patterns {h,v) = (1,2) or {h,v) = (2,1), that is to say 
alternating black and white stripes, are integrable, with quadratic growth of the degrees. 

2x2 patterns. For {h,v) = (2,2) we have different results for the different patterns in 
Figure [3] 



Both [1010] and [0100] are integrable, with quadratic growth of the degrees. 
The calculation of the degrees for [1011] yields the sequence 

{dn} = 1, 2, 4, 8, 18, 41, 93, 215, 493, 1132, 2600, 5970, 13710, 

31487, 72313, 166077, 381417, 875974, 2011788, 4620332 . . . (26) 

This sequence if fitted by the rational generating function 

/ N , n 1 - - 2 + - s6 - + 

g(s) = > dns" = -— -— -, (27) 

^ ^ (l-s)(s + l)(s4-2s3-2s + l)(s4 + l)' ^ ' 

and gives a non vanishing entropy e = log(s) with s the largest root of — 2 s'^ — 2 s + 1, 
approximately e = log(2. 29663). 



Caveat: When computing sequences of degrees, one should in principle consider iterations 
of the whole pattern, but that tends to make the calculations heavier. For the sequence (26), 
this would mean considering only the subsequence formed by odd terms, leading to a growth 
given by the maximal root r of — 4 — 6 — 4 1 + 1 . Of course r = . 



2x3 patterns. We have examined all the period (2, 3) patterns. The various nonequivalent 
patterns are depicted in Figure |4] 

The computations show that what matters is not just the proportion of black and white 
cells, but the actual conformation of the pattern. For example the period (2, 3) patterns 
[010110] and [001011] have an equal number of black and white cells. The first one is 
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Figure 4: Nonequivalent period (2,3) patterns. 

integrable (quadratic growth of the degrees) while the latter is not, as may be seen from the 
sequence: 

{dn} = 1, 1, 2, 3, 5, 9, 19, 41, 84, 169, 329, 631, 1199, 2287, 4412, 8627, 17059, 33941, 

67573, 134071, 264576, 519343, 1015531, 1982461, 3871597, 7574863, 14855790 .. .(28) 

This sequence has exponential growth, but is not long enough to determine an exact value 
of the entropy. The approximate value is log (1.96). 

Out of the 14 period (2, 3) nonequivalent patterns, we have one linear case (all white 
[000000]), eight integrable cases ( [010000], [010100], [110000], [011000], [010101], [010110], 
[111100], [mill]), and four non integrable ones ([110100], [001011], [110101], [111001]). 
The following pictures show the aspect of two integrable cases and two non integrable ones. 





011000 010110 

Figure 5: Two integrable (2, 3) patterns. 



Remark: The entropy calculations for these patterns can be made equally well with the 
relations { white = bottom 1, black=bottom2 } or with { white = topi, black=top2}. Both 
would give the same results. 

From the above results, one may already conclude that random distributions are expected 
to be non-integrable. 
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6 Lax pair for a 2 x 2 sublattice 



We will next push the Lax concept and ZCC to a 2 x 2 sublattice described in Figure [7) 
(Such sub-lattices have been discussed previously, e.g., in [H].) To determine the evolution 
we need 5 initial values, marked with black disks and expect to get values for the vertices at 
the open circles. Since the Lax matrices belong to PGL{2, C) the zero curvature conditions 
can provide at most three equations, and thus if everything works well the evolution is 
determined. It is also clear that this will not work for bigger sub-lattices as we would then 
need to provide more than 3 values. 



M{X22,X2 



M{x2,x) 



L{Xi22, X22) L(^Xii22, 3^122) 

-^(3^11225 3^112) 

M(xii2,Xii) 



• ( 


3 Y 


1) 

L (^Xi , x^ 




L{xii, 3^1)^^ 



Figure 7: The 2x2 configuration. Values at black discs are initial data (x, Xi, Xn, X2, 3:22)5 
and values at open circles (3:112, 3:122, 3:1122) should be determined by the evolution. 

The zero curvature condition for the 2x2 sublattice is given by 

M(a:ii22 ,3:112) M(xii2, Xn) L{xu,Xi) L{xi,x) ~ 

L(X1122, 3:122) -^^(3:122, 3:22) M{X22, X2) M{x2, x). (29) 
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6.1 Flipped Hie 



In the flipped Hie case the Lax matrices were given in Eq. (22) Now using this on a 2 x 2 
sub-lattice we get just one condition, namely 

A(p, xi22, XU22) + A(p, X22, 3:122) + A(g, X2, X22) + A(g, x, X2) = 
A(g, X112, X1122) + A(g, Xn, 0:112) + A(p, Xi, Xu) + A(p, x, Xi). 

From this one can in principle solve X1122 in terms of the other variables. However, this 
equation does not determine the values for X122 or X112 and therefore these Lax matrices fail 
to give the evolution. 



6.2 HI 



For this basic model the Lax matrices were given in Eq. (14). Condition (29) leads to 
equations that have two rational solutions: The regular one 

{p - q){xi - X2) 



X122 



Xi + 
X2 + 



x\\22 = x + (p-q) 
and an exotic solution 



(xi - X2)(x - Xii) -{p-qY 
{q-p){x2 - xi) 

(x2 - xi)(x - X22) - {q-pY 

{p - g)(xii + X22 - 2x) + 2(x - xii)(x - X22)(xi 



X2) 



{p - qY - (X - Xii)(x - X22)(Xl - X2)^ 



X1122 

X122 
X112 



Xll + X22 



X, 



X112 - Xi + X2, 



Xl - 



(xi - X2)[{P - r){x - X22) + (g - r){x - Xn)] 



{p - r)(x - X22) - {(I - r){x - Xll) - (x - Xii)(x - X22)(xi - X2) 



(30a) 
(30b) 
(30c) 



(31a) 
(31b) 

,(31c) 



The regular solution could also be obtained using the evolution on the original lattice, 
first solving for X12. As a consequence X122 depends only on x,xi,X2,X22 and X112 only 
on x,Xi,X2,Xii. The exotic solution is different, as X122 and X112 both depend on all initial 
values. Furthermore, it depends on p — r and q — r, and not solely on p — q, as is the case 
for the regular solution. 

We may view the variables x, Xn, X22, X1122 as associated to the vertices and Xi, X2, X112, X122 
as associated to the bonds of the 2x2 sublattice. 

In the algebraic entropy analysis the vertex variables are linear and for the bond variables 
we find the sequence of degrees 



{dn} = 1, 4, 13, 28, 49, 76, 109, 148, 193, 244,301, 364,433,... 
This sequence can be fitted with the generating function 

1 + 4 + s 



:i-s) 



(32) 



(33) 
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The sequence has quadratic growth, signahng integrabihty. 



What it the nature of the exotic solution? Since the vertex variables have independent 
linear evolution we can solve the equation with X2n,2m = F{n) + G{m). When this is 
substituted into the bond equations they give a non-autonomous generalization of a Yang- 
Baxter map: using the coarse grained indexing Wn,m = X2n,2m, Xn,m = X2n+i,2m, yn,m = 
X2n,2m+i, 1-6., Xi = X, X2 = Y, Xn2 = ^1, 3:122 = -^2 we havc 

Y^-X = P{X,Y), X2-Y = P{X,Y). (34a) 

The solution w = constant is not allowed and if either F or G is constant, P collapses to 
P = ±{x — y). In the generic case, denoting 



f{n) 



p — r 



g[m) 



q — r 



we get 



p {X-Y)\j{n)^g{m)\ 



X-Y-f{n)+g{m)- 

After the further translation X X + f{n) + T{n,m), Y i-> Y + g{n) + T{n,m), where T is 
a solution of T(n, m + 1) — T{n, m) = 2g{m), T{n + 1, m) — T{n, m) = 2f{n) we finally get 



f{n f -g{mf 
X -Y 



(34b) 



which is a non-autonomous version of the Adler map [21J (aka Fy in the classification |22j). 
The situation can be described by the following diagram: 



standard solution ^ 



HI Lax for 2x2 sublattice 



HI 



exotic solution (31) 



non-autonomous Fy (34) 



6.3 H3 

The phenomenon described in the previous section is not generic. Indeed, applying the same 
coarse-graining to an arbitrary integrable quad-equation will lead to a system having only 
one rational solution (the regular one coming from the original lattice). 

We have, however, found more examples where an exotic rational solution exists. Here is 
one, provided by the lattice modified KdV (ImKdV) (aka H35=o)- In that case the defining 
relations of the exotic solution are: 

3^1122 a; = a;iiX22, (35a) 
0:122X1 = Xii2a;2, (35b) 
X122 a;ii2 {q'^X22 - r'^x){x - Xu)pxi - {p^Xu - r'^x){x - X22)qx2 



X2 Xi {r'^X22 — q'^x){xii — x)px2 — (r^xn — p'^x){x22 — x)qxi 



(35c) 
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In the algebraic entropy analysis the sequence of degrees for the vertex variables has linear 
growth as expected, while the sequence for the bonds is the same as for HI (see above). 

Now the equation on the vertex variables (35a) can be solved with 

X2n,2m = F{n)G{m), 

and if we introduce 

, r'^F(n + l) -p'^F(n) , , r'^Gim ^\) - (f'Gim) 



F(n + 1)-F(n) ' ^ ' G(m + 1) - G(m) 

we obtain the bond equations in the form 

X2 _ Fi _ (g^ + - g{m))pX - (/ + - f{n))qY 
Y ~ X ~ fin)qX - g{m)pY ' 



(36) 



using the previously introduced notation. With the further scaling 

X{n, m) (— )■ pT{n, m) X{n, m)/f{n), Y{n, m) 1— )■ qT{n, m) Y{n, m)/ g{m), 
where T solves 

T{n,m + 1) /T{n,m) = {q'^+r'^ — g{m))/g{m), T{n+l,m) /T{n,m) = {p^ + r'^ — f{n)) / f{n), 



equation (36) reduces to 

X.y p, y^.JLp, (37) 

a{nj p[mj X — Y 

where 

a{n) = XpV[f{n){f{n) - p' - r% /3(n) = \q' /[g{n){g{n) - q' - r% 

This is nothing but a non-autonomous version of Fju in the classification of [22]. The 
diagram presented for HI works also for H3((5 = 0). 

We have found that this phenomenon occurs also for the lattice modified KdV and for 
the lattice Schwarzian KdV. We have examined in some detail the properties of the models 
defined in this way and we will present these results elsewhere [16]. 



7 Discussion 



We have discussed the strength of the Lax pairs (or the zero curvature condition) and BT 
through various examples. We have found several cases where the ZCC does not uniquely 
determine the evolution but allows two possibilities. This happens, e.g., in the fiipped Hie 
model. If one then builds an infinite lattice by arbitrarily choosing for each cell one of the two 
allowed relations, the result is sometimes integrable and sometimes shows nonzero entropy. 
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If the zee is pushed to a 2 x 2 sublattice we get more examples where it is ambiguous 
and yields both the regular solution as well as an exotic one. The latter cannot be generated 
by some equation in the sublattice, because in the exotic solution the variables X122 and Xiu 
depend on both xn and 0:22 which is not possible using the equations on the elementary 
squares. 

The equations we have obtained in this way can be interpreted as having vertex and edge 
variables: the variables with even number of indices live at the vertices while the ones with an 
odd number of indices live on the edges. The edge variables evolve as in a non-autonomous 
functional Yang-Baxter equation. For more details see [16]. 

Acknowledgments 

One of us (JH) was partially supported by the Ville de Paris through the program "Research 
in Paris 2010". Some of the computations were done using REDUCE [23j. 

References 

[1] FW Nijhoff, A Ramani, B Grammaticos and Y Ohta, On Discrete Painleve Equations 
Associated with the Lattice KdV Systems and the Painleve VI Equation, Stud. Appl. 
Math. 106 (2001) 261. 

[2] FW. Nijhoff and A. Walker, The discrete and continuous Painleve VI hierarchy and the 
Gamier systems, Glasgow Math. J. 43 A, 109 (2001). 

[3] A.I. Bobenko, Yu.B. Suris. Integrable systems on quad-graphs. Internat. Math. Res. 
Notices, 11, 573-611 (2002). 

[4] V Adler, A Bobenko and Yu Suris, Classification of Integrable Equations on Quad- 
Graphs. The Consistency Approach, Commun. Math. Phys. 233 (2003) 513. 

[5] FW. Nijhoff, Lax pair for the Adler (lattice Krichever-Novikov) system, Phys. Lett. 
A297 49-58 (2002). 

[6] J. Atkinson, J. Hietarinta and FW. Nijhoff, Seed and soliton solutions of Adler 's lattice 
equation, J.Phys.A: Math. Theor. 40 (2007) F1-F8. 

[7] J. Atkinson, J. Hietarinta and FW. Nijhoff, Soliton solutions for Q3, J. Phys. A: Math. 
Theor. 41 (2008) 142001 (11 pp). 

[8] J. Hietarinta and D.J. Zhang, Soliton solutions for ABS lattice equations: II Casoratians 
and bilinearization, J. Phys. A: Math. Theor. 42 (2009) 404006 (30pp). 

[9] V. E. Adler, A. I. Bobenko, and Yu. B. Suris. Discrete Nonlinear Hyperbolic Equations. 
Classification of Integrable Cases, Funct. Anal. App., 43, 3-17 (2009). 



18 



J. Atkinson. Bdcklund transformations for integrable lattice equations J. Phys. A: Math. 
Theor. 41 135202 (2008). 

P D Xenitidis and V G Papageorgiou. Symmetries and integrability of discrete equations 
defined on a black-white lattice, J. Phys. A: Math. Theor. 42 (2009) 454025 (13pp.). 

D. Levi and R. I. Yamilov. On a nonlinear integrable difference equation on the square, 
Ufimsk. Mat. Zh. 1, 101-105 (2009). 



R. Boll. Classification of 3D consistent quad- equations, arXiv : 1009 . 4007v2 [nlin . SI] 
3 Nov 2010. 

F. Calogero and M.C Nucci, Lax pairs galore, J. Math. Phys. 32, 72 (1991). 

Mike Hay, Discrete Lax Pairs, Reductions and Hierarchies, Thesis (University of Sydney, 
2008). 

J. Hietarinta and C. Viallet, Integrable lattices with vertex and bond variables. 

M. Bellon and C-M. Viallet. Algebraic Entropy Comm. Math. Phys. 204 (1999) 425-437. 



C-M. Viallet Algebraic entropy for lattice equations, arXiv:math-ph/0609043 



S. Tremblay, B. Grammaticos, and A. Ramani, Integrable lattice equations and their 
growth properties, Phys. Lett. A 278 (2001) 319-324. 

J. Hietarinta and C. Viallet. Searching for integrable lattice maps using factorization, J. 
Phys. A: Math. Theor. 40 (2007) 12629-12643 

V.E. Adler, Recuttings of polygons, Funct. Anal. App., 27, 141-143 (1993). 

V.E. Adler, A.I. Bobenko and Yu.B. Suris, Geometry of YangBaxter Maps: pencils of 
conies and quadrirational mappings, Commun. Anal. Geom., 12, 967-1007 (2004). 

[23] A. Hearn, REDUCE User's Manual Version 3.8 (2004) 



jhttp : / / reduce-algebra . sourcef orge . net/ 



19 



